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ABSTRACT 

We present the discovery of four absorption lines in the X-ray spectrum of 
the Seyfert Galaxy NGC 1365, at energies between 6.7 and 8.3 keV. The lines 
are detected with high statistical confidence (from > 20a for the strongest to 
~ 4<t for the weakest) in two XMM-Newton observations 60 ksec long. We also 
detect the same lines, with lower signal-to-noise (but still > 2a for each line) 
in two previous shorter (~ 10 ksec) XMM observations. The spectral analysis 
identifies these features as FeXXV and FeXXVI Ka and K/3 lines, outflowing 
with velocities varying between ~ 1000 to ~ 5000 km s _1 among the observations. 
These are the highest quality detections of such lines so far. The high equivalent 
widths (EW(i^a) ~ 100 eV) and the Ka/ 'K/3 ratios imply that the lines are due 
to absorption of the AGN continuum by a highly ionized gas with column density 
Nh ~ 5 x 10 23 cm -2 at a distance of ~ 50-100 Rs from the continuum source. 

Subject headings: Galaxies: AGN — X-rays: galaxies — Galaxies: individual 
(NGC 1365) 

1. Introduction 

Warm absorbers are common among Active Galactic Nuclei (AGNs). About half of 
the X-ray spectra of local bright Seyfert galaxies show evidence of a warm absorber with 

1 Harvard- Smithsonian Center for Astrophysics, 60 Garden St. Cambridge, MA 02138 USA 
2 INAF - Osservatorio di Arcetri, L.go E. Fermi 5, Firenze, Italy 

3 XMM-Newton Science Operations Center, European Space Astronomy Center, ESA, Apartado 50727, 
E-28080 Madrid, Spain 

4 Dipartimento di Fisica, Universita degli Studi "Roma Tre", Via della Vasca Navale 84, 1-00146 Roma, 
Italy 



- 2 - 



column densities in the range 10 22 -10 23 cm -2 (e.g., Reynolds 1997). High resolution grating 
observations performed with Chandra and XMM-Newton have shown that the warm absorber 
has a typical temperature of ~ 10 6 K, and, in some cases, outflowing velocities of a few 
10 3 km s" 1 (Sako et al. 2001, Kaspi et al. 2002, Krongold et al. 2003). Recently more 
extreme, "hot" absorbers, have been discovered through the observation of He-like and H- 
like iron absorption lines in several AGNs (e.g., in NGC 3783, Kaspi et al. 2002, Reeves 
et al. 2004, in PG 1115+080, Chartas et al. 2003, in PG 1211+143, Pounds et al. 2003, in 
MCG-6-30-15, Young et al. 2005). 

Here we report the discovery of four strong absorption lines in each of four XMM- 
Newton observations of the Seyfert Galaxy NGC 1365, in the 6.7-8.3 keV band, which we 
identify as FeXXV and FeXXVI Ka and K(3 lines. NGC 1365 (z=0.0055, de Vaucouleurs 
et al. 1992) is an X-ray absorbed Seyfert Galaxy, with a 2-10 keV flux of the order of 
0.5-2xl0" n erg s" 1 cm" 2 (Risaliti et al. 2000, 2005, hereafter R05), which shows extreme 
X-ray variability: changes from reflection-dominated to transmission dominated states are 
observed on time scale as short as 3 weeks (R05). When in a transmission dominated state, 
the measured cold column density varies between ~ 1 and ~ 5 x 10 23 cm" 2 . 



2. Data Analysis and Results 

Two XMM-Newton (Jansen et al. 2001) observations of 60 ksec each were obtained in 
December 2003 and July 2004. Two shorter (~ 10 ksec) observations were obtained earlier 
in 2003. A spectral analysis of all these observations showed the source to be in transmission 
dominated state (R05). The observation log is reported in Table 1, together with the main 
parameters of the continuum fit which will be described in the next Section. 

We reduced the EPIC PN and MOS data using the SAS 6.0 package. We extracted 
a spectrum from a circular region of 30" radius, and the background from nearby regions 
free from bright serendipitous sources. In all cases the background counts are negligible 
at all energies (on average ~ 2% of the signal). In particular, we checked that no narrow 
background line 1 is present in the energy range which is most crucial for our analysis, i.e., 
between 6.7 and 8.3 keV, where the absorption lines are detected. The data obtained with 
the two MOS were merged. Calibration matrices were computed for each spectrum. We 
rebinned the spectra in order to have at least 20 counts per bin. This allows us to use x 2 
minimization in the model fitting. 



1 e.g., from Ni, Cu, Zn, see http://xmm.vilspa.esa.es/external/xmm_user_support/documentation/uhb/node35.html 
(Lumb et al. 2002). 
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For all four spectra a good fit is obtained with a model consisting of an absorbed 
power law with Nh between 1 and 5 x 10 23 cm~ 2 , plus a cold reflection component (using 
the PEXRAV model in XSPEC, Magdziarz & Zdziarski 1 995) and a narrow emission line at 
6.4 keV. At energies below the photoelectric cut-off a thermal emission component dominates, 
with kT ~ 0.8 keV (R05). A broad emission line, probably relativistic, is also needed to 
properly fit the 4-6 keV region 2 . In addition to these "standard" components, four absorption 
lines in the energy range 6.7-8.3 keV are strongly requested by the fit. A complete analysis 
of the emission spectrum of NGC 1365 will be presented elsewhere, together with a detailed 
study of the spectral variability during the two long observations. Here we concentrate on 
the most striking features in these spectra, i.e., the group of absorption lines present in the 
6.7-8.3 keV spectral region. In order to visually show the relevance of these lines, in Fig. 1 
we plot the 5-10 keV residuals to a model with the same components as above, except for 
the four absorption lines, and fitted ignoring the energy ranges 6.7-7.2 keV and 7.8-8.3 keV. 
The spacing of the lines strongly suggests that they are FeXXV and FeXXVI Ka and Kf3. 

Accordingly, in the global 0.5-10 keV fits the four lines were fitted with narrow Gaus- 
sians, with free energies allowed for each single line and each single fit. All the four lines 
are significantly detected in all spectra, except in two cases where the detection is marginal 
(FeXXV K/3 in OBS 1 and OBS 2). For the two highest quality observations, OBS 3 and 
OBS 4, the fits were performed on PN and MOS data separately. 

Since the best fit energies are close to those of the four lines mentioned above, we 
estimated the blueshift /redshift velocities corresponding to the differences between the mea- 
sured and theoretical energies for these lines. The velocities are shown in Fig. 2. In both 
observations the results are in agreement with a common outflow velocity for the four lines. 
Therefore, we performed a new fit fixing the energies of the four lines to the theoretical 
values (taking into account for the redshift of the source), and allowing for a single free 
redshift /blueshift. We obtained as good a fit as in the previous case where all the lines were 
left free (A(x 2 ) < 8 in the global fit in all observations, with three less free parameters). 

The best fit line parameters are listed in Table 2. The continuum parameters, listed in 
Table. 1, are typical of X-ray spectra of obscured, Compton-thin Seyfert galaxies, and will 
be discussed in detail in a forthcoming paper. 

In Fig. 3 we show the 5-10 keV spectrum and best fit model for the two long observations. 
The measured Ka/Kf3 ratios (Table 3) are significantly smaller than those expected from 



Alternative continuum models are possible, for example with the inclusion of partial covering. Here we 
do not investigate this point since we are only interested in an analytic fit to the continuum, in order to have 
a correct determination of the narrow absorption lines parameters. 
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the oscillator strengths ratios 3 . 

In all four spectra the lines are blueshifted with respect to the expected theoretical 
line energy peaks. The velocity shifts are between 1,000 and 5,200 km s" 1 (Table 2), and 
are statistically highly significant, except for the first observation, where the data are in 
agreement with zero velocity. The differences between the outflow velocities are highly 
significant, and cannot be due to instrumental effects. In order to check the reality of this 
shift, we compared the low energy thermal emission lines in the ~ 1 keV region (which are 
expected to be constant, given the good S/N and the extent of the thermal emission region 
resolved by Chandra, R05) and found that they overlap perfectly (A(E) < 10 eV). 

3. Discussion 

The main results of our spectral analysis are the following: 

• We detected four absorption lines in the energy range 6.7-8.3 keV in four XMM-Newton 
observations of NGC 1365, with equivalent width in the 50-150 eV range. We identified these 
lines as FeXXV and FeXXVI Ka and Kf3. 

• In all observations the lines are blueshifted. Blueshift velocities are the same for all 
lines in the same observation, but vary among the different observations from ~ 1, 000 to 
~ 5,000 km s" 1 . 

• The Ka/K/3 ratios are much smaller than the theoretical oscillator strength values (Ta- 
ble 3), and are consistent with staying constant between the observations. 

In the following we discuss the composition (ionization state, column density), the 
geometrical structure, and origin of the hot absorber. 

Composition. The low Ka/Kfi ratios can only be explained if saturation effects are 
significant. Any other physical explanation is ruled out, since (a) the theoretical ratio is 
well established, since it involves transitions to the fundamental level in relatively simple 
(He-like and H-like) atomic species, and (b) the energy difference between the lines is too 
small for the continuum slope to play a significant role: for a photon index T = 2.0 we 
have Fe(Q-7 keV)/Fe(7.9 keV)= 1.17. Following B05, an equivalent width of the Ka lines 
EW~100-150 eV can be obtained with Nh > 10 23 cm" 2 , a turbulent velocity of the absorbing 
gas v t = 1,000 km s _1 , and an ionization parameter logUx = 0. 4 We note that the EPIC 



3 Throughout this paper we refer to the atomic physics values of the NIST Atomic database 
(http://physics.nist.gov), and reported in Table 1 of Bianchi et al. 2005, hereafter B05. 

4 We adopt the definition in B05: Ux — (J^ j^dv)/{A-Kr 2 cn e ), similar to the one first introduced by 
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resolution at 8 keV corresponds to velocities of ~ 5, 000 km s _1 . With the available statistics 
we find that the minimum detectable velocity is of the order of 2,000 km -1 . Therefore, no 
useful constraint on v t can be obtained from our fits. 

We calculated the curve of growth for the Kf3 lines using the same model as B05. In 
order to have a precise estimate of the EW ratios, we also considered the contribution of 
the FeXXV line at 8.29 keV, which is blended with the FeXXVI K(3 in low resolution 
spectra, and has an oscillator strength about half of that of the FeXXVI Kf3. A further 
minor contribution could come from the NiXXVII (He like) Ka at 7.98 keV. Its oscillator 
strength is 0.19 that of FeXXV Kf3. An overabundance of nickel could make this contribution 
important. However, we note that in this case the line profile would be significantly altered, 
given that the energy separation between the two lines (AE = 80 eV) is comparable to the 
EPIC resolution (~ 150 eV for both PN and MOS cameras; but note that with the high 
statistics available the peak energies of the lines is determined with a much better precision, 
~10— 15 eV). This is not observed (Fig. 1). We conclude that the contribution of nickel lines 
cannot be significant. In Fig. 4 we plot the curve of growth for the FeXXV and FeXXVI Ka 
lines (from B05) and the Ka/Kj3 ratio. Note that both the large EWs and the low Ka/Kf3 
ratios require an absorber with a column density Nh ~ 5 x 10 23 cm~ 2 and a turbulent 
velocity a ~500-l,000 km s" 1 . The right panels of Fig. 4 show the importance of saturation 
effects in the Ka/ Kf3 ratios. High column densities are required in order to reproduce the 
small observed ratios. 

The above estimates are computed adopting a fixed ionization parameter Ux- Therefore 
we expect they slightly change for different choices of Ux- However we do not expect to find 
acceptable solutions for significantly lower or higher ionization parameters: if (a) \ogUx » 
the iron atoms are almost completely stripped, making impossible to obtain as high EWs 
as observed; if (b) —1 < \ogU x < the He-like Fe atoms are overabundant with respect to 
the H-like ones, because the ionizing continuum is enough to strip all the L shell electrons, 
but not to strip the 2 inner electrons. The ratio between FeXXV and FeXXVI lines would 
then be much higher than observed. This argument is discussed in more detail in B05. 

The turbulent velocity is significantly higher than the thermal velocity of a gas with the 
same ionization state (<Jt ~100-200 km s _1 ). As a consequence, the line widths probably 
indicate true bulk motions. These might be related to the observed changes in outflow 
velocity. 

Geometrical properties. The high ionization parameter required to explain the line 



Nctzcr (1996). For ease of comparison, if \ogUx — 0, the ionization parameter denned as U' x = Lx/(n e r 2 ) 
has a value U' x ~ 5, 000 (for an ionizing continuum with T = 2). 
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ratios can be used to put constraints on the distance of the line absorber from the X- 
ray continuum source, assuming that the ionization state is due to UV radiation from the 
central source, rather than to a high thermal temperature of the gas. For logUx = 0, Nh = 
5 x 10 23 cm -2 , and a 2-10 keV photon index T = 2 we obtain R ~ 2 x 10 15 (AR/R) L 42 cm, 
where L42 is the 2-10 keV luminosity in units of 10 42 erg s _1 (Table 1) and AR is the thickness 
of the absorber along the radial direction. R is the average distance of the absorber, so we 
expect AR/R ~ 1. We note that a distance of ~ 10 15 cm is of the order of what a gas with 
a velocity v ~ 4, 000 km s -1 (intermediate between the two measurements in the last two 
observations, OBS 3 and OBS 4) would cover in ~ 6 months between OBS 3 and OBS 4 
(~ 2 x 10 7 sec). This is in qualitative agreement with our assumption of an absorbing 
gas with a thickness of the same order of the distance from the central source. A smaller 
thickness (AR/R < 1) would imply an even smaller distance from the central source. 

Origin. The short distance inferred from the high ionization state of the absorber, and 
the rough black hole mass estimate (Mbh ~ 10 8 M , R05) imply that the gas is located at 
the radii of the accretion disk, ~ 50 (M/Mg) -1 Schwarzschild radii from the black hole. Hot 
gas could be present in this region either as the external part of the hot corona believed to 
be responsible of the X-ray emission in AGNs 5 (e.g., Haardt & Maraschi 1991), or as the 
inner part of a wind arising from the disk (Murray & Chiang 1995, Proga et al. 2000). 

We note that the observed blueshift velocity is smaller than the escape velocity at 
50 {M/M s )- l R S) which is of the order of about 50,000 (M/M 8 )-°- 5 km s _1 . However, the 
observed velocity is of the order of that predicted in the vertical part of the funnel-shaped 
wind of Elvis (2000). The main problem with the wind scenario is that radiation force in 
not an effective acceleration mechanism in our case, since the gas is overionised (so line 
absorption is not an efficient way to transfer momentum from radiation to gas) and the 
bolometric luminosity (of the order of a few 10 43 erg s _1 adopting a standard X-ray to 
bolometric correction, Risaliti & Elvis 2004) is a small fraction of the Eddington luminosity, 
L>edd ~ 2 x 10 46 (M/Mg) _1 erg s _1 , making Thomson scattering acceleration negligible with 
respect to gravity. However, two mechanisms could provide the observed highly ionized gas: 
(a) a poloidal magnetic field could effectively extract the gas from the accretion disk and 
form a wind (Blandford & Payne 1982, Konigl & Kartje 1994); (b) gas arising from the 
accretion disk due to local instabilities would be ionized by the central UV/X-ray source, 
and initially expand towards the vertical direction (with respect to the disk plane) due to its 
internal pressure. The gas would then fail to form a wind if no external force is present, but 
a fraction of the line of sight, at small angles with respect to the disk plane, would be covered 



5 The inner corona, emitting the bulk of X-ray radiation, has to have kT ~ 100 eV, therefore it is too hot 
to absorb its own radiation. 
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by such gas, with the required ionization state and column density, and with a (temporary) 
outflowing velocity of the order of a few 10 3 km s _1 . Similar properties are predicted for the 
inner "hitchhiking gas" in radiation driven wind models (Murray & Chiang 1995). 

We are grateful to F. Nicastro and N. Brickhouse for useful discussions. This work was 
partially supported by NASA grants NAG5-13161, NNG04GF97G, and NAG5-16932. 
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Fig. 1. — Contributions to x 2 i n the 5-10 keV range for the four MOS observations 
of NGC 1365. The continuum model has been fitted ignoring the energy intervals 6.7- 
7.2 keV and 7.8-8.3 keV. The four vertical lines show the rest frame energy of the four lines 
FeXXV Ka, FeXXVI Ka, FeXXV K/3, FeXXVI K(3. 
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Fig. 2. — Best fit velocities for the four iron absorption lines in the two highest quality 
observations, OBS 3 and OBS 4. For each line we separately plot the results from the PN 
spectra (blue circles) and the MOS spectra (red squares). The vertical solid line shows the 
best fit value obtained requiring a common velocity for the four lines with ±1 a values shown 
as vertical dashed lines. 




Fig. 3. — Hard X-ray (5-10 keV) spectra and best fit model for the two long XMM-Newton 
observations of NGC 1365. The background spectrum is shown for each observation. 
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Fig. 4. — Left panels: curves of growth for the FeXXV and FeXXVI Ka lines, for log Ux = 0. 
Right panels: Ka/K/3 ratios for FeXXV and FeXXVI. The yellow shaded regions are the 
90% confidence intervals for the averages obtained from the four observations. 
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Table 1: Observation Log and Spectral Parameters 



OBS 


Date 


t a 


Count s b 


F c 


L d 


r 


N H e 


1 

2 
3 
4 


2003 Jan 16 

2003 Aug 13 

2004 Jan 17 
2004 Jul 24 


14.8 
7.1 
59 
60 


14866 
9466 
121967 
74056 


0.48 
0.73 
1.37 
0.76 


1.0 
1.4 
1.7 
0.8 


o i cr+0.03 
z - id -0.06 
9 09+O.I2 
z - oz -0.02 
9 cr+0.05 
z - od -0.03 
9 1 1 +0.03 
z - iJ --0.04 


51+1 
37 + 2 

i6.3i8;| 

29.418;| 



a Net observing time in ksec. 

^Source counts in the 0.5-10 keV interval. 

c 2-10 keV observed flux, in units of 10~ n erg s _1 cm -1 . 

^Intrinsic 2-10 keV luminosity in units of 10 42 erg s _1 . 

e In units of 10 22 cm -2 . 
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Table 2: NGC 1365 - Data on absorption lines 



Line 


E(keV) 


OBS 1 


OBS 2 


OBS 3 


OBS 4 






EW (eV)' 


a A(x 2 )C 


EW (eV) a 


A(x 2 ) c 


EW (eV) a 


A(x 2 ) c 


EW (eV) a A(x 2 ) c 


FeXXV Ka 


6.697 


nu_ 70 


17 


135±30 


30 


120±20 


173 


154±20 264 


FeXXVI Ka 


6.966 


04+50 
y ^-50 


11 


111±30 


20 


ml:™ 


52 


1451^ 167 


FeXXV Kf3 


7.880 


81±55 


5 


100±60 


6 


93±10 


43 


108±20 60 


FeXXVI K(3 d 


8.268 


127±60 


12 


130±60 


10 


80±15 


18 


110±25 59 


V(km s^ 1 ) 13 




iooo: 


hlOOO 
-1500 


20001 


1800 
400 


25001 


500 
300 


c-900+400 
OZUU_ 400 



"Errors are at a 90% confidence level. 
b Outflowing velocity of the absorbing gas. 

c A(x 2 ) with respect to the best fit model without the four absorption lines. The four lines are forced to have 
the same blueshift, so the model including the lines has 5 more parameters than the no-line model. 
d The FeXXVI K(3 line is blended with the FeXXV Kj (see text for details). 
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Table 3: FeXXV and FeXXVI Ka/Kp ratios 





OBS 1 


OBS 2 


OBS 3 


OBS 4 


pa 


Th. 6 


FeXXV 
FeXXVI 


1.4±1.7 
0.7±0.7 


1.4±1.1 
0.9±0.6 


1.4±0.4 
1.4±0.7 


1.4±0.4 
1.3±0.4 


1.4±0.3 
1.1±0.3 


5.0 
3.2 



a : Average ratios: b : Oscillator strengths ratio. In R2 the contribution of FeXXV K7 has been taken into 
account. 



